Oscillations and defect turbulence in a shallow fluidized bed 
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We report an experimental study of the dynamics of an air-fluidized thin granular layer. Near- 
onset behavior of this shallow fluidized bed was described in the earlier paper (Tsimring et al, 
1999). Above the threshold of fluidization the system exhibits a Hopf bifurcation as the layer starts 
to oscillate at a certain frequency due to a feedback between the layer dilation and the airflow 
drag force. After application of temporal band-pass filtering of this frequency we discovered the 
spatio-temporal dynamics in the form of defect turbulence. This type of dynamics is natural for 
spatio-temporal systems close to the threshold of a Hopf bifurcation. At high flow rates, low- 
frequency short-wavelength structures appear in addition to the long-wavelength excitations. A 
simple model describing the instability and occurrence of oscillations in a shallow fluidized bed, is 
proposed. 
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The non-equilibrium dynamics of granular materials 
have been a subject of growing attention among physi- 
cists in the last few years Granular materials re- 
veal a host of interesting phenomena (compaction, self- 
organized criticality, patterns, convection, etc.) when are 
subjected to external driving. A typical way of driving 
the granular medium out of equilibrium is to fluidize it 
with an external gas or liquid flow. Such "fluidized beds" 
are widely used in the industry for mixing solid and liq- 
uid chemicals One of the biggest problems associ- 
ated with this use of fluidized beds is an instability which 
leads to macroscopic inhomogeneities and fluctuations in 
the bulk (bubbling, slugging). At large flow rates, this 
leads to a developed three-dimensional turbulent regime 
resembling boiling liquid. 

In the recent paper ||, we proposed to study the dy- 
namics of air-fluidized granular materials in a bed with 
a large aspect ratio between the horizontal size and the 
thickness of the layer. In this geometry, the dynamics 
become quasi-two-dimensional and can be investigated 
by imaging the surface of the granular layer. In Ref. [||, 
we studied the onset of fluidization and found that it has 
the features of a phase transition, with very long tran- 
sient fluctuations just below the fluidization threshold. 
In this paper, we report the results of the experimental 
study of the thin granular layer dynamics at high airflow 
rates, above the threshold of fluidization. We find that 
above this threshold the layer begins to oscillate at a cer- 
tain frequency which depends on the layer thickness and 
the airflow rate. Processing of surface images revealed 
the occurrence of propagating waves which form small 
disordered short-lived spirals. At larger airflow rate, ir- 
regular small-scale oscillating cellular pattern appears on 
the surface of the bed. The temporal frequency of these 
oscillating bubbles is significantly less than that of large- 
scale waves. We believe that these small-scale pertur- 
bations are generated by the large amplitude, long- wave 



layer oscillations in a way similar to the mechanism of 
Faraday instability in fluid or granular layers subject to 
external vertical vibrations (see, e.,g., (|J^]). The mech- 
anism of layer oscillations is related to the dependence 
of the drag force from the air flow acting on the granu- 
lar layer, on the volume fraction of particles in the layer. 
When the layer is dilated, the drag force diminishes, and 
the layer falls back due to gravity. This picture leads us 
to the formulation of the simple dynamical model which 
qualitatively agrees with our experimental results. 

The experimental setup is similar to the one used in 
Ref. but it features a larger round porous bronze 
plate (diameter 15 cm, thickness 0.5 cm, average pore size 
6/im). As a fluidizing fluid we use compressed dry air. 
The granular material consists of monodisperse spherical 
bronze particles of size 0.15 mm, and the bed depth was 
varied from 3 to 7 particles deep. The patterns on the 
surface of the granular layer were illuminated using the 
low-angle light from a circular fluorescent light around 
the bed, and recorded using high-speed Kodak SR-C dig- 
ital CCD camera. We also used a pressure transducer 
SenSym SX01DP1 to record the airflow fluctuations near 
the free surface of the granular layer. 

When the flow rate is below certain threshold value, 
the layer remains completely static. The value of this 
critical airflow depends sensitively on the grain parame- 
ters (size, shape, density) as well as on the thickness of 
the layer. Unless explicitly noted otherwise, we will be re- 
ferring to the bed thickness 0.65 mm. For that thickness, 
the critical airflow speed was 12.1 cm/sec. Slightly above 
the threshold, the fluidization of the layer is highly non- 
uniform. As described in our earlier paper Ref. ||, the 
fluidization occurs in small localized regions surrounded 
by completely static grains. 

At higher values of the airflow the whole layer becomes 
fluidized and exhibits noticable vertical oscillations. The 
power spectrum of the pressure fluctuations near the 
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free surface of the layer is shown in Figure |T|. We mea- 
sured the peak frequency and the magnitude as a func- 
tion of the airflow velocity (Figure ||). The amplitude 
of the oscillations increases linearly and the frequency / 
slightly decreases as the airflow speed increases. Figure 
||a shows dependence of the peak frequency on the air- 
flow speed for different values of layer thickness. A naive 
scaling of variables by the gravity g and the thickness 
of the layer ho (/(ho/g) 1 ^ 2 , V(gho)~ 1 ^ 2 ) does not con- 
verge these data to a single curve. It indicates that addi- 
tional non-dimensional parameters involving the particle 
size d and the air viscosity v play an important role in 
the mechanism of these oscillations. For example, if one 
chooses scaling Vdj (h^g) 1 / 2 for the airflow, then the data 
lie rather close to a single curve (see Fig.||b). 
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FIG. 1. Power spectrum of the pressure fluctuations near 
the layer at the airflow rate 15.8 cm/sec. The arrow indicates 
the peak corresponding to the large-scale layer oscillations. 
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FIG. 2. Amplitude (a) and frequency (b) of the spectral 
peak in the power spectrum of layer oscillations as a function 
of the airflow rate V for the layer thickness 0.65 mm. Linear 
regression of the amplitude (dashed line in a) yields the os- 
cillation threshold of 11.2 cm/s, which is slightly lower than 
the fluidization threshold 12.1 cm/s 
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FIG. 3. a - frequency / of the spectral peak in the power 
spectrum of layer oscillations as a function of the airflow 
rate V; b - rescaled frequency f(hfg) 1 ^ 2 as a function of 
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FIG. 4. Spatial structure of the 0.65 mm layer oscillations: 
a - snapshot of the granular layer surface for V =15.8 cm/sec; 
b - corresponding snapshot from the temporal sequence fil- 
tered by the narrow-band filter with central frequency / = 35 
Hz; c - spatial distribution of the phase of 35 Hz oscillations; 
d - snapshot from the temporal sequence for V =42.2 cm/sec 
filtered by a low-pass filter with the cutoff frequency 25 Hz. 
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High-speed imaging of the surface at 500 fps also con- 
firms the occurrence of the layer oscillations. In order 
to determine the spatial structure of these oscillations, 
we performed temporal filtering of the recorded image 
sequences. We used a band-pass filter to extract os- 
cillations at the peak frequency / for each pixel, and 
then combined these filtered signals back to create the 
sequence of temporally-filtered images ||. 

A raw snapshot from the video sequence taken at the 
airflow speed 15.8 cm/sec, and the corresponding filtered 
snapshot are shown in Figure ^a,b. Viewing filtered im- 
ages in sequence [^| reveals that the oscillations form 
propagating waves, moreover, these waves are organized 
in small irregular spirals. The sense of direction of the 
wave propagation can be determined from the spatial dis- 
tribution of the phase of oscillations at the filtered fre- 
quency / (Figure |]c). Singularities of the phase corre- 
spond to the cores of the spirals, or topological defects. 
The phase distribution changes over time, as the defects 
move, annihilate, and new defects are born. Such "defect 
turbulence" is to be expected near the onset of the Hopf 
bifurcation in spatially extended systems. It is generi- 
cally described by the Ginzburg-Landau equation for the 
complex amplitude of the oscillation |lo|. Similar de- 
fect turbulence has been observed in spatially-extended 
Belousov-Zhabotinsky reaction jnj]. 

At higher air flow, in addition to large-scale layer os- 
cillations, irregular small-scale perturbations of the layer 
appear. These perturbations correspond to the low- 
frequency part of the layer oscillations spectrum. Figure 
||d shows the snapshot of the surface image at the airflow 
42.2 cm/sec obtained by low-pass filtering of the image 
sequence with the cut-off frequency 25 Hz. The spatial 
power spectrum of these perturbations averaged over an- 
gle is shown in Figure || together with the corresponding 
spectrum for the peak frequency 31 Hz. As one can see, 
spatial spectrum of low-frequency oscillations exhibit a 
broad peak at wavenumber k ~ 1.6 cm -1 , whereas the 
spectrum of 31 Hz oscillations decays with the wavenum- 
ber. 
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FIG. 5. Azimuthally- averaged spatial power spectrum of 
oscillations at the dominant frequency 31.2 Hz (dashed line) 
and low-frequency oscillations / < 25 Hz (solid line) for the 
airflow speed 42.2 cm/sec 




The mechanism of the shallow fluidized bed oscilla- 
tions can be described as follows. The layer dynamics 
are determined by the interplay of two forces: gravity 
and the drag force from the gas flow. While gravity g 
remains constant, the drag force varies depending on the 
volume fraction of particles in the layer. When the layer 
lies on the surface and is densely packed, the drag force 
is maximal. On the other hand, when the layer is lifted 
in the air and dilated, it becomes more "transparent" for 
the airflow, and the drag force decreases. This effect is 



known for deep fluidized beds |12|. In the simplest form, 
we could assume that the drag force T is inversely pro- 
portional to the height of the layer ft, T = 7/ ft, where 
constant 7 depends on the nominal flow speed V, close- 
packed thickness of the layer fto, gas viscosity, and the 
properties of the individual particles (size, shape). Then 
the second Newton's law yields 
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where subscript t denotes time derivative, and p is the 
total mass of grains per unit area. In non-dimensional 
variables H = ft/fto and T = (g/fto) 1 / 2 ^, this equation 
reads 



Hi 1 = — - 1, 



(2) 



where P = j/pgho- The amplitude and frequency of 
periodic oscillations for this nonlinear oscillator equation 
depend on the initial conditions. Assuming that the layer 
always comes down to the surface where it has a minimal 
thickness fto, we impose initial conditions H(0) = 1 and 
Ht(0) = 0. With these initial conditions, amplitude A 
and frequency Q = uj(hofg) 1 ^ 2 of oscillations are uniquely 
determined by the parameter P. For P slightly above 1, 
the amplitude of oscillations A w P — 1 and the frequency 
is ft w P -1 / 2 . So, in qualitative agreement with the 
experiment, near the threshold the amplitude increases 
approximately linearly with P, and the frequency slowly 
decreases with P. 

While the model Eq.(|l|) explains the nature of the layer 
oscillations, it does not address the mechanism of the in- 
stability which leads to the oscillations. Indeed, since 
there is always some dissipation due to friction and in- 
elastic particle collisions, the oscillations would decay, 
and the stable equilibrium bed thickness at which the 
drag force balances gravity, would be reached. To over- 
come dissipation, an additional instability mechanism is 
needed to excite the layer oscillations. The mechanism of 
this instability is related to the known mechanism of in- 
stability for deep fluidized beds |l^| , in which it is caused 
by the inertia of particles in the gas flow and leads to 
the time delay between the volume fraction variations 
and the drag force changes. Similarly, we can assume 
that the drag force T is not enslaved to the thickness the 
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layer ft, but approaches 7/ft asymptotically. This addi- 
tional degree of freedom can be described by the equation 
aT t = 7ft -1 — T, where parameter a controls the rate of 
inertia in the dependence between T and ft. We shall also 
modify the equation for the layer thickness dynamics (jl]) 
to include a thickness-dependent dissipation. Indeed, the 
amount of dissipation varies strongly as a function of the 
layer thickness, as it affects the frequency of inelastic 
collision among grains. As the thickness approaches the 
close-packed limit, the number of collisions, and the dis- 
sipation rate, diverge. We assume that this effect is cap- 
tured by an additional dissipative term —crht/(h — ho) 2 
(the actual power of singularity at ft = fto in not impor- 
tant here). Now, the extended model in non-dimensional 
form reads 
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where F = T/gp,a! = a(g/h ) 1/2 , a' = a(^ftg)- 1 / 2 . Ad- 
ditional relaxation dynamics described by Eq. ([|) directly 
leads to an instability of the thin fluidized layer. Lineariz- 
ing Eqs.(||),(||) near the fixed point F = 1, H = P yields 
the dispersion equation for perturbations ex exp(— zOT), 
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where we dropped primes at the non-dimensional param- 
eters a and a. At small a and a, the complex frequency 
of layer oscillations is given by 
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As one can see from this formula, the instability always 
occurs at large enough P. Notice that for small a and 
cr, the solution of Eqs.(^),(||) at large T approaches that 
of Eq.(|) with initial conditions if (0) = 1; H T (0) = 0. 
This is not surprising, since these conditions are enforced 
by diverging dissipation rate at H = 1, and during the 
rest of the oscillation cycle, non-conservative corrections 
are negligibly small. 

To conclude, in this paper we studied the dynamics 
of a shallow fluidized granular layer. Experiments with 
air-driven shallow fluidized bed showed that above fru- 
idization threshold, the layer exhibits oscillations at a 
certain frequency. These oscillations are in fact propa- 
gating waves which have a structure of small disordered 
spirals. At larger airflow speeds, in addition to the high- 
frequency spiral waves, low- frequency short- wavelength 
perturbations are observed in the bed. It is feasible that 
these small "bubbles" are formed via subharmonic exci- 
tation by the primary high-frequency oscillations, as in 
the Faraday instability. We proposed a simple dynami- 
cal model describing high-frequency oscillations in which 



the layer is driven by the airflow drag force and gravity. 
The drag force in turn depends on the volume fraction of 
particles in the layer, and therefore on the height of the 
layer surface. The underlying assumptions for this model 
are based on experimental facts and the intuitive phys- 
ical picture, and a systematic derivation of the model 
from the first principles may yield more complicated re- 
lations among the layer parameters. Still, this model 
describes the instability leading to the layer oscillations, 
and agrees on a qualitative level with the observed de- 
pendence of the oscillation magnitude and frequency on 
the airflow rate. A more detailed study of the model 
including its systematic derivation and generalization to- 
wards spatiotemporal dynamics will be the subject of a 
separate publication. 
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